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7.1 INTRODUCTION 


Flows at high Reynolds number reveal that the viscous effects are confined within 
the boundary layers. Far away from the solid surface, the flow is nearly inviscid 
and in many cases it is incompressible. We now aim at developing techniques for 
analyses of inviscid incompressible flows. 

Incompressible flow is a constant density flow, and we assume p to be 
constant. We visualize a fluid element of defined mass moving along a streamline 
in an incompressible flow. Because the density is constant, we can write 

V-V = 0 (7.1) 

Over and above, if the fluid element does not rotate as it moves along the 
streamline, or to be precise, if its motion is translational (and deformation with no 
rotation) only, the flow is termed as irrotational flow. It has already been shown 
in Sec. 3.3.5 that the motion of a fluid element can in general have translation, 
deformation and rotation. The rate of rotation of the fluid element can be 
measured by the average rate of rotation of two perpendicular line segments. The 
average rate of rotation co, about z-axis is expressed in terms of the gradients of 
velocity components (refer to Chapter 3) as 

dv du 
dx dv 

Similarly, the other two components of rotation are 



dw 

dv 

and 

1 

du 

dw 



Cft, = - 



dy 

dz 


■' 2 

dz 

dx 
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As such, they are components of co which is given by 



In a two-dimensional flow, (0 : is the only non-trivial component of the rate of 
rotation. Imagine a pathline of a fluid particle shown in Fig. 7.1. Rate of spin of 
the particle is co : . The flow in which this spin is zero throughout is known as 
irrotational flow. A generalized statement is more appropriate: For irrotational 
flows, V x V = 0 in the flow field. 



Fig. 7.1 Pathline of a fluid particle 


Therefore for an irrotational flow, the velocity V can be expressed as the gradient 
of a scalar function called the velocity potential, denoted by 0 


V = V0 

Combination of Eqs (7.1) and (7.2) yields 


(7.2) 


V 2 0 = 0 


(7.3) 


From Eq. (7.3) we see that an inviscid, incompressible, irrotational flow is 
governed by Laplace’s equation. 

Laplace’s equation is linear, hence any number of particular solutions of 
Eq. (7.3) added together will yield another solution. This concept forms the 
building-block of the solution of inviscid, incompressible, irrotational flows. A 
complicated flow pattern for an inviscid, incompressible, irrotational flow can be 
synthesized by adding together a number of elementary flows which are also 
inviscid, incompressible and irrotational. 

The analysis of Laplace’s Eq. (7.3) and finding out the potential functions are 
known as potential flow theory and the inviscid, incompressible, irrotational 
flow is often called as potential flow. However, the following elementary flows 
can constitute several complex potential-flow problems 

1. Uniform flow 

2. Source or sink 

3. Vortex 
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7.2 ELEMENTARY FLOWS IN A 
TWO-DIMENSIONAL PLANE 


7.2.1 Uniform Flow 


In this flow, velocity is uniform along y-axis and there exists only one component 
of velocity which is in the x direction. Magnitude of the velocity is U 0 . 

From Eq. (7.2) we can write 


iU 0 + jO 



- d(f) 
J dy 


or 


d<p 

dx 


u 0 


d<p 

~dy =0 


whence <p=U 0 x+C l (7.4) 

Recall from Sec. 4.2.2 that in a two dimensional flow field, flow can also be 
described by stream function y. In the case of uniform flow 

dy dy 

— = f/ 0 and - — 
dy dx 


= 0 


so that l F=U 0 y + K l (7.5) 

The constants of integration Cj and K x in Eqs (7.4) and (7.5) are arbitrary. The 
values of y and (p for different streamlines and velocity potential lines may change 
but flow pattern is unaltered. The constants of integration may be omitted and it is 
possible to write 

y=U 0 y, (/)=U 0 x (7.6) 

These are plotted in Fig. 7.2(a) and consist of a rectangular mesh of straight 
streamlines and orthogonal straight potential-lines. It is conventional to put 
arrows on the streamlines showing the direction of flow. 

In terms of polar (r - 6) coordinate, Eq. (7.6) becomes 

y = U 0 r sin 0, </> = U 0 r cos 0 (7.7) 

If we consider a uniform stream at an angle a to the x-axis as shown in Fig. 7.2b, 
we require that 


u 


U 0 cos a = 


dy 

dv 


d<p 

dx 


and v — Uq sin cc = — = ~~~ (7.8) 

dx dy 

Integrating, we obtain for a uniform velocity U 0 at an angle cqthestream function 
and velocity potential respectively as 
y = U 0 (y cos a - x since), 0 = U 0 (x cos ce + y since) 


(7.9) 
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(a) 



(b) 


Fig. 7.2 (a) Flownet for a uniform stream 

(b) Flownet for uniform stream with an angle a with x-axis 

7.2.2 Source or Sink 

Consider a flow with straight streamlines emerging from a point, where the 
velocity along each streamline varies inversely with distance from the point, as 
shown in Fig. 7.3. Only the radial component of velocity is non-trivial ( v g = 0, 
v,_ = 0). 



Fig. 7.3 Flownet for a source flow 
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Such a flow is called source flow. In a steady source flow the amount of fluid 
crossing any given cylindrical surface of radius r and unit length is constant ( m ) 
77? =2 nr v r p 


or 


v.. = 


1 


2 np r 

where, K is the source strength 

„ _ 777 A 

A — 


_A_ J_ 
In r 


K_ 

r 


(7.10a) 


(7.10b) 


2 np 2 n 

and A is the volume flow rate 

Again recall from Sec. 4.2.2 that the definition of stream function in cylindrical 
polar coordinate states that 

dw 

Ve ~~lh 

Now for the source flow, it can be said that 
l dijf _ K 
/- 39 r 


1 d\jf 
Vr r dO 


and 


(7.11) 


(7.12) 


and 


dip 

dr 


= 0 


Combining Eqs (7.12) and (7.13), we get 
V r = K0+C l 

However, this flow is also irrotational and we can write 




~ c dtb 
i v e = 

dr 


1 d(f) 
~d0 


J 


or 

or 


v r 


dtp 

dr 


and 


Vg= 0 = 


1 dtp 
r dO 


d<p _ _ K 

dr Vr r 


or <j> = K In r + C 2 


(7.13) 

(7.14) 


(7.15) 


Likewise in uniform flow, the integration constants C 1 and C 2 in Eqs (7.14) and 
(7.15) have no effect on the basic structure of velocity and pressure in the flow. 
The equations for streamlines and velocity potential lines for source flow become 
V f = K0 and <p = K\nr (7.16) 

where K is defined as the source strength and is proportional to A which is the rate 
of volume flow from the source per unit depth perpendicular to the page as shown 
in Fig. 7.3. If A is negative, we have sink flow, where the flow is in the opposite 
direction of the source flow. In Fig. 7.3, the point O is the origin of the radial 
streamlines. We visualize that point O is a point source or sink that induces radial 
flow in the neighbourhood. The point source or sink is a point of singularity in the 
flow field (because v r becomes infinite). It can also be visualized that point O in 
Fig. 7.3 is simply a point formed by the intersection of plane of the paper and a 
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line perpendicular to the paper. The line perpendicular to the paper is a line 
source, with volume flow rate (A) per unit length. However, for sink, the stream 
function and velocity potential function are 

t j/ = -K6 and 0 = K In r (7.17) 

7.2.3 Vortex Flow 

In this flow all the streamlines are concentric circles about a given point where 
the velocity along each streamline is inversely proportional to the distance from 
the centre, as shown in Fig. 7.4. Such a flow is called vortex (free vortex) flow. 
This flow is necessarily irrotational. 


Y=- 



2k 


In ?•] 


In a purely circulatory (free vortex flow) motion, we can write the tangential 
velocity as 

Circulation constant 
v e = - 


r/27r 

v e= - 

r 

where T is circulation, 

Also, for purely circulatory motion one can write 
v r = 0 

With the definition of stream function, it is evident that 

d Y „ - 1 d V 


(7.18) 


(7.19) 


Vn = ~ 


and 


v = 


dr ' r d9 

Combining Eqs (7.18) and (7.19) with the above said relations for stream 
function, it is possible to write 

y/= In r+C l 
2 n 


(7.20) 
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Because of irrotationality, it should satisfy 


iv r + jv 


e 


' d(j) ~ 1 d(j) 

dr ^r d6 


Eqs (7.18) and (7.19) and the above solution of Laplace’s equation yields 

<p = JL e + C 2 (7.21) 

2 n 

The integration constants C, and C 2 have no effect whatsoever on the structure of 
velocities or pressures in the flow. Therefore like other elementary flows, we 
shall consistently ignore such constants. It is clear that the streamlines for vortex 
flow are circles while the potential lines are radial. These are given by 

w = —— In r and 6 (7.22) 

2 n 2 n 

In Fig. 7.4, point O can be imagined as a point vortex that induces the circulatory 
flow around it. The point vortex is a singularity in the flow field ( v e becomes 
infinite). It is also discerned that the point O in Fig. 7.4 is simply a point formed 
by the intersection of the plane of a paper and a line perpendicular to the plane. 
This line is called vortex filament of strength T, where T is the circulation around 
the vortex filament and the circulation is defined as 


r= <j)Vds 


(7.23) 


In Eq. (7.23), the line integral of the velocity component tangent to a curve of 
elemental length d.v, is taken around a closed curve. It may be stated that the 
circulation for a closed path in an irrotational flow field is zero. However, the 
circulation for a given path in an irrotational flow containing a finite number of 
singular points is constant. In general this circulation constant Y denotes the 
algebraic strength of the vortex filament contained within the closed curve. 
From Eq. (7.23) we can write 

r = cj)V ds = f (udx + vdy + wdz) 

For a two-dimensional flow 


T= j> (udx + v dy) 

or r = <j) Ceos a ds (7.24) 

Consider a fluid element as shown in Fig. 7.5. Circulation is positive in the anti¬ 
clockwise direction (not a mandatory but general convention). 


, fa c 
U+ dy Sy 

v I Sy 

I 

u 

Fig. 7.5 Circulation in a flow field 











The McGraw-Hill Companies 


284 


Introduction to Fluid Mechanics and Fluid Machines 


or 


or 


8T = u 8x + 

sr = 


v- 


dv „ 
v —Sx 
dx 

Sx Sy 


Sy- 


du „ ^„ 

-dv ox-vov 

dy ' J 


dv du 
dx dy 
8T = 2co z 8A 

or 8T/8A = 2co : = Q, (7.25) 

Physically, circulation per unit area is the vorticity of the flow. 

Now, for a free vortex flow, the tangential velocity is given by Eq. (7.18) as 

r/ 27 r c 

v g = -= — 

r r 

For a circular path (refer Fig. 7.5) 

a = 0 , 


T7- C 

V=v e = — 


2 K 


Thus, 


C 


T= A-rd9 = 2nC 


(7.26) 


It may be noted that although free vortex is basically an irrotational motion, 
the circulation for a given path containing a singular point (including the origin) 
is constant (2 n C ) and independent of the radius of a circular streamline. 
However, if the circulation is calculated in a free vortex flow along any closed 
contour excluding the singular point (the origin), it should be zero. Let us look at 
Fig. 7.6 (a) and take a closed contour ABCD in order to find out circulation about 
the point, P around ABCD 


1 ABCD - V e AB ' 1 

There is no radial flow 


r i - Vr BC 0 2 - c) + v 6rn r 2 dfl + v rnA (r 2 - ;q) 


r DA 


c 


v,. = v,. = 0, v n = — and v ft = 

r BC , oa ’ Vjr - Vrn 


c 


ABCD 


8cd 


-c 


c 


= -qd0 H— —r^dO — 0 


(7.27) 


If there exists a solid body rotation at constant (0 induced by some external 
mechanism, the flow should be called a forced vortex motion (Fig. 7.6b) and we 
can write 

v e = cor and 


2 K 

r= <§>v e ds = ^(Or rdO — 2nr 2 (0 (7.28) 

o 

Equation (7.28) predicts that the circulation is zero at the origin and it 
increases with increasing radius. The variation is parabolic. 

It may be mentioned that the free vortex (irrotational) flow at the origin 
(Fig. 7.6a) is impossible because of mathematical singularity. However, 
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physically there should exist a rotational (forced vortex) core which is shown by 
the dotted line. Below are given two statements which are related to Kelvin’s 
circulation theorem (stated in 1869) and Cauchy’stheoremonirrotational motion 
(stated in 1815) respectively 

(i) The circulation around any closed contour is invariant with time in an 
inviscid fluid. 

(ii) A body of inviscid fluid in irrotational motion continues to move 
irrotationally. 



(b) 


Fig. 7.6 (a) Free vortex flow (b) Forced vortex flow 
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7.3 SUPERPOSITION OF ELEMENTARY FLOWS 

We can now form different flow patterns by superimposing the velocity potential 
and stream functions of the elementary flows stated above. 

7.3.1 Doublet 

In order to develop a doublet, imagine a source and a sink of equal strength K at 
equal distance s from the origin along x-axis as shown in Fig. 7.7. 

From any point P{x, y) in the field, r, and r 2 are drawn to the source and the 
sink. The polar coordinates of this point (r, 0) have been shown. 



Fig. 7.7 Superposition of a source and a sink 

The potential functions of the two flows may be superimposed to describe the 
potential for the combined flow at P as 

<t> = K\nr l -K\nr 2 (7.29) 

Similarly 

y/ = K(6 l - 0 2 ) = -Ka (7.30) 

where, a = (0 2 - 0j) 

We can also write 


tan 6\ = —— 


x + s 


y 

and tan 0 2 = —‘— 


(7.31) 


= -Jr 2 + s 2 + 2rs cos 0 


and 




r 2 + s 2 - 2 rs cos 0 


Now using the above mentioned relations we find 
tan (0 2 - 0,) = 


or 


tan a = 


tan0 2 — tan0[ 

1 + tan0 2 tan0j 

yx + ys — _vx + ys 

2 2 
X — S 


1 / 

( 2 A 

/ 


J/ 

l x ~s ) 


(7.32) 
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or 


tan a = 


2ys 


2,2 2 

x + y — s 


(7.33) 


Hence the stream function and the velocity potential function are formed by 
combining Eqs (7.30) and (7.33), as well as Eqs (7.29) and (7.32) respectively 

f \ 


y/ = -K tan 1 


2 vs 


2.2 2 
\x +y -s 


K, 


' r 2 + s 2 + IrscosO ^ 
y r 2 + s 2 — IrscosOj 


(7.34) 

(7.35) 


Doublet is a special case when a source as well as a sink are brought together in 
such a way that s —> 0 and at the same time the strength K ( A/2n) is increased to 
an infinite value. These are assumed to be accomplished in a manner which makes 

the product of s and — (in limiting case) a finite value X. Under the aforesaid 
n 

circumstances 


y- 


2 vs 


A 

2n x 2 + y 2 - s 2 


[Since in the limiting case tan 1 a = a\ 

V -ysin0 

V=-X- , ' , - 

x~ + y r 

From Eq. (7.35), we get 


(7.36) 


(j) = [In (r 2 + s 2 + 2rs cos 0) - In (r 2 + s 2 - 2rs cos 0)] 
An 


, A 

or <p= — 
An 


In | [r 2 + s 1 j 


f 


2rscos0 


r 


■ ln<^ (r 2 + s 2 ) 


2rscos0 

2.2 
r + s 


, A 
or m= — 
An 


2rscos 1 
r 2 + s 2 2 


^2rscos 0 

2 ~, 2 ~ 
\ r +S J 


2rscos 0 

2 ~, 2 ~ 

r +s 




2 


3 


2 rs cos 0 1 

2rscos6 

1 

2rscos 0 



r 2 +s 2 2 

2 . 2 

V r ) 

3 

2 , 2 
r +s 




, A 

or 0= — 
An 


4rscos0 2 
3 


f 2rs cos 0^ 


2 , 2 
r +5 


2 . 2 
r + s 


In the limiting condition the above expression can be written as 

%r cos 6 


0 : 


2 , 2 
r +5 
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or 


</>: 


£COS0 


(7.37) 


We can see that the streamlines associated with the doublet are 
Xsmd _ ^ 

L i 

r 

If we replace sin0 by y/r, and the minus sign be absorbed in C h we get 


x 4 = q 

r 

In terms of cartesian coordinate, it is possible to write 


2,2 
x +y - 


■-%-y = 0 

Q 


(7.38a) 


(7.38b) 


Equation (7.38b) represents a family of circles. Forx= 0, there are two values 
of y, one of which is zero. The centres of the circles fall on the y-axis. On the 
circle, where y = 0, x has to be zero for all the values of the constant. It is obvious 
that the family of circles formed due to different values of C, must be tangent to 
x-axis at the origin. These streamlines are illustrated in Fig. 7.8. Due to the initial 
positions of the source and the sink in the development of the doublet, it is certain 
that the flow will emerge in the negative x direction from the origin and it will 
converge via the positive x direction of the origin. 



Fig. 7.8 Streamlines and velocity potential lines for a doublet 













The McGraw-Hill Companies 


Flow of Ideal Fluids 


0*9 


However, the velocity potential lines are 


^COS0 

r 


= Ki 


In cartesian coordinate this equation becomes 


x 2+y!_X x= 0 (7.39) 

A 1 

Once again we shall obtain a family of circles. The centres will fall on x-axis. For 
y = 0 there are two values ofx, one of which is zero. Whenx = 0, y has to be zero 
for all values of the constant. Therefore these circles are tangent to y-axis at the 
origin. The orthogonality of constant (// and constant (j) lines are maintained as we 
iron out the procedure of drawing constant value lines (Fig. 7.8). In addition to 
the determination of the stream function and velocity potential, it is observed 
from Eq. (7.37) that for a doublet 


dd) 

Vr ~lfr 


-ycosO 


(7.40) 


As the centre of the doublet is approached, the radial velocity tends to be infinite. 
It shows that the doublet flow has a singularity. Since the circulation about a 
singular point of a source or a sink is zero for any strength, it is obvious that the 
circulation about the singular point in a doublet flow must be zero. It follows that 
for all paths in a doublet flow F = 0 

r=cfv-ds=0 (7.41) 

Applying Stokes Theorem between the line integral and the area-integral 

JJ,V x V) cL4 = 0 (7.42) 

From Eq. (7.42), the obvious conclusion is V x V = 0, i.e., doublet flow is an 
irrotational flow. 

At large distances from a doublet, the flow approximates the disturbances of a 
two dimensional airfoil. The influence of an airfoil as felt at distant walls may be 
approximated mathematically by a combination of doublets with varying 
strengths. Thus the cruise conditions of a two dimensional airfoil can be simu¬ 
lated by the superposition of a uniform flow and a doublet sheet of varying 
strengths. 


7.3.2 Flow About a Cylinder Without Circulation 

Inviscid-incompressible flow about a cylinder in uniform flow is equivalent to the 
superposition of a uniform flow and a doublet. The doublet has its axis of 
development parallel to the direction of the uniform flow. The combined potential 
of this flow is given by 


r 


(7.43) 
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and consequently the stream function becomes 


v= u dy- 


ysin 0 


(7.44) 


r 

In our analysis, we shall draw streamlines in the flow field. Intwo-dimensional 
flow, a streamline may be interpreted as the edge of a surface on which the 
velocity vector should always be tangent and there is no flow in the direction 
normal to it. The latter is identically the characteristics of a solid impervious 
boundary. Hence, a streamline may also be considered as the contour of an 
impervious two-dimensional body. Figure 7.9 shows a set of streamlines. The 
streamline C-D may be considered as the edge of a two-dimensional body while 
the remaining streamlines form the flow about the boundary. 



Fig. 7.9 Surface streamline 


Now we follow the essential steps involving the superposition of elementary 
flows in order to form a flow about the body of interest. A streamline has to be 
determined which encloses an area whose shape is of practical importance in 
fluid flow. This streamline will describe the boundary of a two-dimensional solid 
body. The remaining streamlines outside this solid region will constitute the flow 
about this body. 

Let us look for the streamline whose value is zero. Thus we obtain 


c oV _Z 

r 

replacing y by r sin 0, we have 


sin 0 | u o r % 
r J 


= 0 


(7.45) 


(7.46) 


If 9 = 0 or 0 = n, the equation is satisfied. This indicates that the x-axis is a part 
of the streamline y/= 0. When the quantity in the parentheses is zero, the equation 
is identically satisfied. Hence it follows that 


C V/2 

X_ 

\ U oJ 


It can be said that there is a circle of radius 


/ V /2 

X 


(7.47) 


which is an intrinsic part of 


the streamline y/ = 0. This is shown in Fig. 7.10. Let us look at the points of 
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intersection of the circle and x-axis, i.e. 
the points A and B. The polar coordinates 
of these points are 


r = 


r = 


f V/2 

X_ 

\ U oJ 

( V /2 

X_ 

\ u oJ 


9= n, for point A 


0=0, for point B 



Fig. 7.10 Streamline \|/ = 0 in a 
superimposed flow of 
doublet and uniform stream 


The velocity at these points are found 
out by taking partial derivatives of the 
velocity potential in two orthogonal 
directions and then substituting the proper values of the coordinates. Thus 

dtp _ tt Q ;^cos0 
dr 


v„ = 


= U 0 cos 0 


(7.48a) 


1 dtp . ^sin0 

y e = = — U 0 sin 0-— 

r dO r 


(7.48b) 


At point A 


At point B 



r \ 

1/2' 

0 = n, r = 

X 






v r = 0,v e =0 



( \ 

1/2" 

0 = 0, r — 

X 






v r = 0, v e = 0 

The points A and B are clearly the stagnation points through which the flow 
divides and subsequently reunites forming a zone of circular bluff body. 

The circular region, enclosed by part of the streamline y/=0 could be imagined 
as a solid cylinder in an inviscid flow. At a large distance from the cylinder the 
flow is moving uniformly in a cross-flow configuration. 

Figure 7.11 shows the streamlines of the flow. The streamlines outside the 
circle describe the flow pattern of the inviscid irrotational flow across a cylinder. 
However, the streamlines inside the circle maybe disregarded since this region is 
considered as a solid obstacle. 



Fig. 7.11 Inviscid flow past a cylinder 
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7.3.3 Lift and Drag for Flow Past a Cylinder Without 
Circulation 


Lift and drag are the forces per unit length on the cylinder in the directions normal 
and parallel respectively, to the direction of uniform flow. 

Pressure for the combined doublet and uniform flow becomes uniform at large 
distances from the cylinder where the influence of doublet is indeed small. Let us 
imagine the pressure p (} is known as well as uniform velocity U (j . Now we can 
apply Bernoulli’s equation between infinity and the points on the boundary of the 
cylinder. Neglecting the variation of potential energy-between the aforesaid point 
at infinity and any point on the surface of the cylinder, we can write 


Po + Up _ Pb + Ub 

Pg 2 g pg 2 g 


(7.49) 


where, the subscript b indicates the surface on the cylinder. As we know, since 
fluid cannot penetrate the solid boundary, the velocity U h should be only in the 
transverse direction, or in other words, only v e component of velocity is present 
on the streamline \jf= 0. 


Thus at 


/ V /2 

X 

XoJ 


U h = 


v e 


at r=(x!U„f 2 


ld± 

r d9 


at r=(xlU„f 2 


(7.50) 


= - 2 Ur. sin 6 


From Eqs (7.49) and (7.50) we obtain 


Pb= Pg 


Ul 

2 g 


Po_ 

Pg 


(2U 0 sin 6Y 

2g 


(7.51) 


The drag is calculated by integrating the force components arising out of pressure, 
in the x direction on the boundary. Referring to Fig. 7.12, the drag force can be 
written as 


U 0 


y 



Fig 7.12 Calculation of drag on a cylinder 
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J I^g3 


D = 


2 K 

J Pb cos 9 


( v /2 

X_ 

K U o ) 


d9 


or 


D — 


2k f \ 

1/2 - 

f l x 


J PS\ 77" 


o WoJ 



Up | Po ( 2U 0 smO ) 2 
2 g PS 


2n 


Dm - 


Similarly, the lift force 


+ pUo a 4 sin 2 

2 


2 n 

f \ 

- [ p h sin 0 

X 

0 

l^o J 


2g 


( \l/2 

x_ 

V^oy 


dO 


cos 0 d0 


cos0 dO 


(7.52) 


(7.53) 


The Eqs (7.52) and (7.53) produced = 0 and L= 0 after the integration is carried 
out. 

However, in reality, the cylinder will always experience some drag force. This 
contradiction between the inviscid flow result and the experiment is usually 
known as D ’Almbertparadox. The reason for the discrepancy lies in completely 
ignoring the viscous effects throughout the flow field. Effect of the thin region 
adjacent to the solid boundary is of paramount importance in determining drag 
force. However, the lift may often be predicted by thepresenttechnique. We shall 
appreciate this fact in a subsequent section. 


7.3.4 Flow About a Rotating Cylinder 

In addition to superimposed uniform flow and a doublet, a vortex is thrown at the 
doublet centre. This will simulate a rotating cylinder in uniform stream. We shall 
see that the pressure distribution will result in a force, a component of which will 
culminate in lift force. The phenomenon of generation of lift by a rotating object 
placed in a stream is known as Magnus effect. The velocity potential and stream 
functions for the combination of doublet, vortex and uniform flow are 

y QOS 0 r 

0 = U 0 x+ - 0 (clockwise rotation) (7.54) 

r 2 n 

'YSVCiQ J -1 

y/ = U 0 y -1-In r (clockwise rotation)(7.55) 

r In 

By making use of either the stream function or velocity potential function, the 
velocity components are 

^■ = 

v e = 




1 dp 

r dO 


U + — 
U 0 + 2 


sin 0 


r 

2nr 


(7.56) 

(7.57) 
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Implicit in the above derivation are x = r cos 6 and y = r sin 9. At the stagnation 
points the velocity components must vanish. From Eq. (7.56), we get 

X 


cos 


e\u, 


'0 2 I — 0 

r 


From Eq. (7.58) it is evident that a zero radial velocity component may occur 

( \ 1/2 

X 


n 


at 9 = ± — and along the circle, r = 
transverse velocity requires 


V^oy 


(7.58) 

/ occur 

. Eq. (7.57) depicts that a zero 


sin 9 = 


—T/2nr 

U 0 +{x!f 2 ) 


or 9 = sin 


-T!2nr 


U n 


2L 

2 


(7.59) 


However, at the stagnation point, both radial and transverse velocity components 
must be zero. 

So, the location of stragnation point occurs at 


r = 


( \ xn 
\ U oJ 


and 


or 


9 = sin 


9 = sin 


f / v /2>i 

X 


2 n 


V^o j 


U 0 + Xj 

-r 


f x_^ 

\ u oj 


2n 


or 


9 = sin 


1 x} 

yU 0J 


-r 


1/2 


2 C/„ 


\l/2 


(7.60) 


4 n(XU 0 ) 

There will be two stagnation points since there are two angles for a given sine 
except for shf 1 (±1). 

The streamline passing through these points may be determined by evaluating 
yrat these points. Substitution of the stagnation coordinate (r, 9) into the stream 
function (Eq. 7.55) yields 


y- 


U n 


\ U 0J 


" X " 

1/2 

X 

sin sin 1 

-r 

+ r in 

f X \ 


f X ' 

1/2 

4n(xU 0 ) 112 

2 n 

y U 0y 


, 1/2 
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v- 


{U 0 X) 112 -(U 0 X) 


1/2 


4 n(XU 0 ) 


1/2 


r 

+ — 

In 

f *1 

2n 




1/2 


or |/ ln 


‘ x A 


V^o j 


1/2 


(7.61) 


Equating the general expression for stream function to the above constant, we 


get 


U 0 r sin 9 - 


_£sin0 T 


A -In r = — In 

r 2 n 2 n 


1 X A 


1/2 


By rearranging we can write 
sin 9 


U 0 r~ — 


2 n 


In/- ± In 


f V /2 

x_ 

V^o j 


\ U oJ 


= 0 


(7.62) 


All points along the circle r = 


f V /2 

X_ 

\ u oj 


satisfy Eq. (7.62), since for this value of 


r, each quantity within parentheses in the equation is zero. Considering the 
interior of the circle (on which y/= 0) to be a solid cylinder, the outer streamline 
pattern is shown in Fig. 7.13. 



A further look into Eq. (7.60) explains that at the stagnation point 


9= sin 


or 


9= sin 


-l 


—(I727t) 

Wo) 1/2 

-(I727T) 

2 U 0 r 


The limiting case arises for ^ = 2, where 9 = sin 1 (-1) = 

U 0 r 

stagnation points meet at the bottom as shown in Fig. 7.14. 


(7.63) 
90° and two 
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Stagnation point 


Fig. 7.14 Flow past a circular cylinder with circulation value 


Y/2n 

U 0 r 


= 2 


However, in all these cases the effects of the vortex and doublet become 
negligibly small as one moves a large distance from the cylinder. The flow is 
assumed to be uniform at infinity. We have already seen that the change in 
strength T of the vortex changes the flow pattern, particularly the position of the 
stagnation points but the radius of the cylinder remains unchanged. 

7.3.5 Lift and Drag for Flow About a Rotating Cylinder 

The pressure at large distances from the cylinder is uniform and given by p 0 . 
Deploying Bernoulli ’ s equation between the points at infinity and on the boundary 
of the cylinder, 


Pb = PS 


2 g 


PO¬ 

PS 


K 

2 g 


(7.64) 


The velocity U b is as such 

_ 1 dtp 


U b = = -2£/ o sin0- — 

r dO In 


Hence, 

From Eqs (7.64) and (7.65) we can write 


Ui > 
X 


1/2 


(7.65) 


Pb = PS 


Ul + Po_ 
2 g pg 


-2U 0 sin0 - 


2n 


Uo 
{ X 


2 g 


(7.66) 


The lift may calculated as (refer Fig. 7.12) 


2k 


L = “I 


Pb 


sin 0 


U n 


d 0 
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O2Z 


or 


2k 


L = - 


put 


+ Po 


-i2 


-2 U n sin 6 ■ 


r u^ 


2 K 


or 


X 

U n 


(sin 0) d 6 


2k 


L = - 


put 


( \ 

X 


\Uq j 


sin 6 + p 0 


( \ 

X 


\Uq j 


sin 14 ut sin 2 0 


or 


+ - 


4 C/ n r sin 9 ( U 0 A 


2 n 

L = -S 

o 


2k 


pu~ Q 


1 

2 r 2 

+ ■ 


( \ 

x 


V X ) 
1 


4 n~ 


\U 0 j 


sin 9 + p 0 


X 


( \ 

X 


( x A 


KU,j 


sin 9 


d 9 


\U 0 j 


sin 9 - 2p U\ sin 3 9 


( 


\ 


U. 


o J 


1 

2 


pu 0 r 

K 


sin 2 9 - 


pr 2 

8 K 2 


V X 2 


sin 9 


d 9 


or L = pU 0 T (7.67) 

The drag force, which includes the multiplication by cos 9 (and integration over 
2n) is zero. 

Thus the inviscid flow also demonstrates lift. It can be seen that the lift 
becomes a simple formula involving only the density of the medium, free stream 
velocity and circulation. In addition, it can also be shown that in two dimensional 
incompressible steady flow about a boundary of any shape, the lift is always a 
product of these three quantities. 

The validity of Eq. (7.67) for any two-dimensional incompressible steady 
potential flow around a body of any shape, not necessarily a circular cylinder, is 
known as the Kutta-Joukowski theorem named after the German fluid dynamist 
Wilhelm Kutta (1867-1944) and Russian mathematician Nikolai J. Joukowski 
(1847-1921). A very popular example of the lift force acting on a rotating body is 
observed in the game of soccer. If a player imparts rotation on the ball while 
shooting it, instead of following the usual trajectory, the ball will swerve in the 
air and puzzle the goalkeeper. The swerve in the air can be controlled by varying 
the strength of circulation, i.e., the amount of rotation. In 1924, a man named 
Flettner had a ship built in Germany which possessed two rotating cylinders to 
generate thrust normal to wind blowing past the ship. The Flettner design did not 
gain any popularity but it is of considerable scientific interest (shown in 
Fig. 7.15). 
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Wind (Vm/s) 



Fig. 7.15 Schematic diagram of the plan view of Flettner's ship 

7.4 AEROFOIL THEORY 

Aerofoils are streamline shaped wings which are used in airplanes and 
turbomachinery. These shapes are such that the drag force is a very small fraction 
of the lift. The following nomenclatures are used for defining an aerofoil (refer to 
Fig. 7.16). 

The chord (c) is the distance between the leading edge and trailing edge. The 
length of an aerofoil, normal to the cross-section (i.e., normal to the plane of a 
paper) is called the span of aerofoil. The camber line represents the mean profile 
of the aerofoil. Some important geometrical parameters for an aerofoil are the 
ratio of maximum thickness to chord {tic) and the ratio of maximum camber to 
chord (h/c). When these ratios are small, an aerofoil can be considered to be thin. 
For the analysis of flow, a thin aerofoil is represented by its camber. 



The theory of thick cambered aerofoils is an advanced topic. Basically it uses 
a complex-variable mapping which transforms the inviscid flow across a rotating 
cylinder into the flow about an aerofoil shape with circulation. 

7.4.1 Flow Around a Thin Aerofoil 

Thin aerofoil theory is based upon the superposition of uniform flow at infinity 
and a continuous distribution of clockwise free vortex on the camber line having 
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circulation density y(s) per unit length. The circulation density 7 ( 5 ) should be 
such that the resultant flow is tangent to the camber line at every point. 

Since the slope of the camber line is assumed to be small, y(.v)ds = y(ri)dri 
(refer Fig. 7.17). The total circulation around the profile is given by 


c 



(7.68) 


0 


A vortical motion of strength ydri at x = )] develops a velocity at the point P 
which may be expressed as 

, y(? 7 ) dn 

dv = - acting upwards 

2n{rj - x ) 

The total induced velocity in the upward direction at P due to the entire vortex 
distribution along the camber line is 

\/ 


y 



pjp'fey) ^ r(iJ ) ds 


P 0, 0 ) 



Camber line 



U, 


x 



c 



[/«, sin a 


U x cos a 


Fig. 7.17 Flow around thin aerofoil 



(7.69) 
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For a small camber (having small a), this expression is identically valid for the 
induced velocity at P' due to the vortex sheet of variable strength 7 ( 5 ) on the 
camber line. The resultant velocity due to [/«, and v(x) must be tangential to the 
camber line so that the slope of a camber line may be expressed as 

dv LC sin a + v v 

— =-= tan a -\ - 

dr U„ cos a LL cos a 

d y j) 

or — = a -1 -[since a is very small] (7.70) 

dx C/ M 


From Eqs (7.69) and (7.70) we can write 

A. I (I)-*) 

Let us consider an element dx on the camber line. Consider a small rectangle 
(drawn with dotted line) around dx. The upper and lower sides of the rectangle are 
very close to each other and these are parallel to the camber line. The other two 
sides are normal to the camber line. The circulation along the rectangle is 
measured in clockwise direction as 


V 1 dx V 2 dx = ydx [normal component of velocity at the camber 
line should be zero] 

or Vy - V 2 = 7 (7.72) 

If the mean velocity in the tangential direction at the camber line is given by 
V s = (V { + V 2 )/2, it can be rewritten as 


v i = v s + J and V 2 ~V*-\ 

In the event, it can be said that if v is very small \v «t/J, V s becomes equal to 
U m . The difference in velocity across the camber line brought about by the vortex 
sheet of variable strength y(x) causes pressure difference and generates lift force. 


7.4.2 Generation of Vortices Around a Wing 

The lift around an aerofoil is generated following Kutta-Joukowski theorem. Lift 
is a product of p, U r „ and the circulation L. Mechanism of induction of circu¬ 
lation is to be understood clearly. 

When the motion of a wing starts from rest, vortices are formed at the trailing 
edge (refer Fig. 7.18). 

At the start, there is a velocity discontinuity at the trailing edge. This is 
eventual because near the trailing edge, the velocity at the bottom surface is 


* For a given aerofoil, the left hand side term of the integral Eq. (7.71) is a known 
function. Finding out 7 ( 77 ) from it is a formidable task. This exercise is not being 
discussed in this text. Interested readers may refer to the books by Glauert [1] and 
Batchelor [2], If 7 ( 77 ) is determined, the circulation T and consequently the lift L = 
p CL, T can easily be calculated. 
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higher than that at the top surface. This discrepancy in velocity culminates in the 
formation of vortices at the trailing edge. Figure 7.18 (a) depicts the formation of 
starting vortex by impulsively moving aerofoil. However, the starting vortices 
induce a counter circulation as shown in Figure 7.18 (b). The circulation around a 
path ( ABCD ) enclosing the wing and just shed (starting) vortex must be zero. 
Here we refer to Kelvin’s theorem once again. 



Fig 7.18 Vortices generated when an aerofoil just begins to move 


Initially the flow starts with the zero circulation around the closed path. 
Thereafter, due to the change in angle of attack or flow velocity, if afresh starting 
vortex is shed, the circulation around the wing will adjust itself so that a net zero 
vorticity is set around the closed path. 

The discussions in the previous section were for two-dimensional, infinite span 
wings. But real wings have finite span or finite aspect ratio X, defined as 


X= — 
A. 


(7.73) 


where b is the span length and A s is the plan form area as seen from the top. For a 
wing of finite span, the end conditions affect both the lift and the drag. In the 
leading edge region, pressure at the bottom surface of a wing is higher than that at 
the top surface. The longitudinal vortices are generated at the edges of finite wing 
owing to pressure differences between the bottom surface directly facing the flow 
and the top surface (refer Fig. 7.19). This is very prominent for small aspect ratio 
delta wings which are used in high-performance aircrafts as shown in Fig. 7.20. 
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Leading edge 



Fig. 7.20 Counter rotating leading edge vortices generated by a delta-wing 

However, circulation around a wing gives rise to bound vortices that move 
along with the wing. In 1918, Prandtl successfully modelled such flows by 
replacing the wing with a lifting line. The bound vortices around this lifting line, 
the starting vortices and the longitudinal vortices formed at the edges, constitute a 
closed vortex ring as shown in Fig. 7.19. 


Summary 


This chapter has given a brief description of inviscid, incompressible, irrotational 
flows. 

• Irrotationality leads to the condition V x V = 0 which demands V = V0, 
where (f) is known as a potential function. For a potential flow V 2 0 = 0. 

• The stream function i// also obeys the Laplace’s equation V 2 y = 0 for the 
potential flows. Laplace’s equation is linear, hence any number of particular 
solutions of Laplace’s equation added togetherwill yield anothersolution. So 
a complicated flow for an inviscid, incompressible, irrotational condition can 
be synthesized by adding together a number of elementary flows which are 
also inviscid, incompressible and irrotational. This is called the method of 
superposition. 

• Some inviscid flow configurations of practical importance are solved by 
using the method of superposition. The circulation in a flow field is defined 
as T = J V • ds . Subsequently, the vorticity may be defined as circulation per 
unit area. The circulation for a closed path in an irrotational flow field is 
zero. However, the circulation for a given closed path in an irrotational flow 
containing a finite number of singular points is a non-zero constant. 

• The lift around an immersed body is generated when the flow field possesses 
circulation. The lift around a body of any shape is given by L = pU 0 T, where 
p is the density and U 0 is the velocity in the streamwise direction. 


Solved Examples 


Example 7.1 The velocity components of two dimensional incompressible flow 
are u = 2xy andu = a 2 + x 2 - y 2 . Show that a velocity potential function exists and find out 
the velocity potential. 
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Solution The velocity potential function exists only for irrotational flow. The condition 
to be satisfied is 

dv _ du 
dx dv 

Evaluating the derivatives mentioned above, we get 

du 


dv 

-=r— = lx 
dx 


and 


dy 


= 2x 


The flow is irrotational. From definition we can also write 

dip 


also, 


LI -y 

dx 

dip ~ 

dip 

v= IT 
dy 


or <p = xv +f\ (y) 


dip 2,2 2 

or ^ = a + x — y or 

Since both the solutions are same, we can write 


2 2 

0 = a y + x y - 


+/ 2 W 


x 2 y +/i(v) = a 2 v + X 2 y - Y +/ 2 (x) 

,3 

In order to keep the above expression valid for all the values of y, f 2 (x) has to be a 
constant. 


or 


A(y)= ° 2 y- V +/ 2 W 


Thus 


2 2 V 

0 = a v + x y - ^ + constant 


Since 0 = constant and represents a family of lines, 0 may be written without a constant 


0 = a 2 v + jc 2 y - 


Exam ple 7. 2 The flow of an incompressible fluid is defined by u = 2, v = 8x. Does 
a stream function exist? If so, find its expression. 

Solution Compliance of continuity describes the existence of a stream function 
du dv _ d{2) d 


dx dy dx dy 
So, the stream function exists. 

Now we can write 


(8x) = 0 


, dy/ , dy/ , 

d V /= ^ dx + ^7 d - v 

dy/= - vdx + udy 
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or d y/ = - 8xdx + 2d v 

or y/=-4x 2 + 2_v + C 

Dropping the constant C, yt= - 4x 2 + 2 y. 


Example 7.3 Does a velocity potential function (/> = 2(x 2 + 2y - y 2 ) describe the 
possible flow of an incompressible fluid? If so, find out the equation for the velocity 
vector V . Also determine the equation for streamlines. 

Solution For the given 0, in order to describe an incompressible flow, we check with 
the Laplace’s equation 

|f + ||-2(2) + 2(-2)^0 

dx dy~ 

So, a flow field exists. 

The velocity components are 

u = = 2(2x) = 4x 

v= = 2(2-2y) = 4-4y 

Velocity vector V = 4x i + (4 - 4v) j 
Stream function y/ can be expressed as 

dx// ^4x dx + ^y dy 
or dy/= - vdx + udy 

or dy/ = - (4 - 4v) dx + 4x dv 

y/= - J (4- 4y)dx + J 4x dv + C 

y/= - 4x + 4xy + 4 xy + C 

Dropping the constant C, stream function becomes 
y/= 4 (2 xy — x) 


Ex amp le 7.4 The radial velocity of a flow is described by v r 
If v g = 0 at 0=0, find out v e and the stream function for the flow. 


Solution 


\ dy/ k r, 

V r = - =i77 = —j= COS 0 

r r dd 


or 

or 

Now 

We know, 



y/= k -Jr sin 0=/(r) 


v e 


dyt _ k 

dr 247- 


sin 0 + f'(f) 


v g = 0 at 0=0, which depicts 
f'(r) = 0 and / (r) = constant 


-4= cos 0. 
•Jr 
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Therefore 


v g = — 7 = sin 6 and 
2 •Jr 

y/= k -Jr sin 0 


Example 7.5 _ A two dimensional source of volume flow rate A = 2.5 m 2 /s is located 
in a uniform flow ( U 0 ) of 2 m/s. Determine the stagnation point and the maximum 
thickness of resulting half body. 

Solution We have already constructed different flow patterns by superimposing 
elementary flows. An interesting body shape appears if we superimpose a uniform flow 
over an isolated source or sink which is known as Rankine half body (refer Fig. 7.21). 
Let the source be located at the origin. 



(i) Then the stream function of combination is 


¥ =U 0 y + ± tan-'gl 

or 

Y= U 0 v + K tan -1 

Velocity 

..- d W TT , V X 

ll - Uq + K 2 2 

°y x + y 

Similarly 

„ _ _ , v y 

V Jy 2 2 

ax x +y 


At the stagnation point u = 0, V = 0, demand from the above equation y — 0 and 


x = 


A 

2nUg 


2.5 

2nx2 


= -0.2m 


The coordinates of stagnation points are (— a, 0) or (— 0.2, 0) 
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The value of stream function at the stagnation point is 
Y (- 0 . 2 , 0 ) = 0 + tan -1 0 


2 5 , 

or V4ta g = 0 + 2k u = 1 25 m /s 


The half-body is described by dividing streamline 



ii 

n 

= kK 

2K 

or 

u Q y 

+ tan -1 - = 4- 

2k x 2 

or 

u Q y 

A 9 A 

2k ~ 2 or - v 2 U 0 

at 

O 

II 

<£> 

y max = 2 XT ’ max i mum ordinate 

at 

qs 

II 

to|a 

y = , the upper ordinate at the origin 

at 

9= K, 

y = 0 , the stagnation point 

at 

II 

Ha 

y = - 4 , the lower ordinate at the origin 


(ii) However, the equation of the half body becomes 

C/ a v+^tan- 1 ^] = A 

The maximum thickness occurs as x —» °° 

2v+— tan^O = 1.25 

K 

1.25 

y — —=0.625 m 

The maximum thickness = 2y max = 1.25 m 


Ex amp le 7.6 A source at the origin and a uniform flow at 5m/s are superimposed. 
The half-body which is formed has a maximum width of 2 m. Calculate (i) location of 

stagnation point (ii) width of the body at the origin and (iii) velocity at a point |o.7, yj. 


Solution (i) We have seen in Example 7.5, that 


at 

or 

for stagnation point, 


® Tmax 2U 0 2 1 m 

A=2x5xl = 10 nr/s 



j=0 


A 

2nU 0 


-10 
2k x 5 


= -0.32 m 


and 
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(ii) at 

e= 

K 

2 ’ 

y = [from Example 7.5] 

at 

e= 

K 

2 ’ 

y = Xx5 ~ 

0.5 m 

The width of the body 

at the origin 

is 2 x 0.5 = 

- 1 m 

(iii) In polar coordinate 





¥= 

U 0 r sin 6+ K6 

where K = - 


v r = 

1 d¥ _ 

r dd 

= U 0 cos 6 + 

A 

2ky 


v e = 

1 

>s ^ 
ll 

- - U 0 sin 6 


at the point (0.7, k!2) 






v,.= 

A 

2 nr 

10 

2k x 0.7 

2.27 m/s 


v g = - U Q sin 0 = -5 sin y = -5 m/s 
Resultant = ^21? + (5 ) 2 = 5.49 m/s 


Ex amp le 7.7 A line source discharging a flow at 0.6 m 2 /s per unit length is located 
at (^1,0) and a sink of volume flow rate 1.2 m 2 /s is located at (2,0). For a dynamic 
pressure of 10 N/nr at the origin, determine the velocity and dynamic pressure at (1,1). 



P(x,y) 


Fig. 7.22 Source and sink pair 


Solution 


or 


I// at P may be expressed as 
y/=K l 6 l - K 2 6 2 


..._ _A]_ a _ Al a 


2k 


2 


y/= 


0.6 —j 

~2k tan 


x +1 


1.2 . -i 

2/F tan 


dy/ 

_ 0.6 

(x + \) 

1.2 

(x-2) 

dy 

2k 

(x +l ) 2 + y 2 

2k 

(x-2 ) 2 +y 2 


or 
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dy/ 

_ 0.6 

y 

1.2 

y 

dx 

2k 

(x + l)^ + 

2k 

(x-2 f+y 2 


at the origin (0,0) 


0.6 
w = -~— 
2k 

v = 0 


1.2 /_ 2 \ _ 0A 0.6 _ 0.6 
27T \ 4 / 2k + 2k ~ K 


= 0.1909 m/s 


Dynamic pressure 


\pV 2 = 10N/m 2 


or p= =548.3 kg/m 3 

At point (1,1) 

_ 0.6 2 12 ]__(L6 0A 

u 2 k5 + 2k2~5k + 2k 

= 0.0381 + 0.0954 = 0.1335 m/s 

_ 0.6 l 1.2 l _ 0.6 0.6 

L ’ 2k'5 2k'2~\0k + 2k 

= 0.019-0.0954 =-0.0764 m/s 

Resultant = 0.1538 m/s 

Dynamic pressure at (1,1) = j x 548.3 x (0.1538) 2 
= 6.48 N/m 2 


Ex amp le 7.8 _ The wind velocity at a location 5 km away from the centre of a tornado 
(consider inviscid, irrotational vortex motion) was measured as 30 km/hr and the 
barometric pressure was 750 mm of Hg. Calculate the wind velocity 1 km from the 
tornado centre and its barometric pressure. [Density of air = 1.2 kg/m 3 , density of mercury 
= 13.6x 10 3 kg/m 3 ]. 
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Solution p a = 750 mm of Hg = 0.75 m of Hg 


= 0.75 x 13.6x 10 3 x 9.81 = 100.062 kN/m 2 


From free vortex consideration, we can write 

r a V 0a = r b V 6b = C 


at 


£AAA 30x1000 , 

r a = 5000 m, v 6a = 6Q x 6Q = 8.33 m/s 

C= circulation constant = 41650 m 2 /s 
r b = 1000 m, 

v eh = =41.65 m/s 


Bernoulli’s equation between points a and b 

Pa V Qa _ Pb | ^Sb 

pg + 2g ~ pg 2g 

100062 , (8.33) 2 = Pb_ + { 41.65) 2 


or 

or 

or 


1.2x9.81 2x9.81 pg 2x9.81 

— = 8500 + 3.536-88.416 
Pg 

= 8415.12 
Pg 

Pb = 99062 N/m 2 = 99.062 kN/m 2 


Ex amp le 7.9 _ A source with volume flow rate 0.2 m 2 /s and a vortex with strength 
1 m 2 /s are located at the origin. Determine the equations for velocity potential and stream 
function. What should be the resultant velocity at x = 0.9 m and y =0.8 m? 

Solution 

For the source y/= 0 = K x In r 

For the vortex y/ = - K 2 In r, <f> = K 2 d 


Combined 

Combined 

Now 

at 


V'=I¥ 0 -2¥ lnr= i 

^ ln '- + 2¥ 0 4 


v„ = 


In 
_ <90 


0 .10-iln r 
0.1 In r + ^6 


1 


dr 10 nr 

= \d± = ^_ 
6 r dd 2 nr 


x = 0.9 m and y = 0.8 m 
r= <[( 0.9) 2 +(0.8) 2 = 1.204 m 


1 


v r (0-9, 0.8) = 

(0-9, 0.8) = 2nr~2xnx 1.204 


10 xnx 1.204 
1 1 


= 0.026 m/s 
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= 0.132 m/s 

K resultanl = V(0.026) 2 +(0.132) 2 =0.134 m/s 


Example 7.10 A 300 mm diameter circular cylinder is rotated about its axis in a 
stream of water having a uniform velocity of 5 m/s. Estimate the rotational speed when 
both the stagnation points coincide. Estimate the lift force experienced by the cylinder 
under such condition, p of water may be assumed to be 1000 kg/m 3 . 


Solution Stagnation point is given by 

(wTT^) 

When both the stagnation points coincide, the two angles are equal and 0 = -90°. 
Stagnation point is at the lower surface [Fig. 7.14], 


Thus 


An r U, 


= 1 


o 


or T = Anr U Q 

If the cylinder is rotating at an angular speed ft), the circulation due to the equivalent 
forced vortex is 


T= | (wr)rdd=2nr 2 ft) 
2ncor= A nrU 0 
2 U 0 


ft) = 


or 

and 

Lift force 


ft)= = 66.67 rad/s 

T = 47rx0.15x5 
= 9.42 m 2 /s 
=L=pUY 


or 


L = 1000x5x9.42 


L= 47100 N/m 


kg 


As such, the lift is calculated per metre length of the cylinder 


So 


Lift = 47.1 kN/nr 
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Exercises 


7.1 Choose the correct answer for the following questions 

(i) Stream function is defined for 

(a) all 3-D flow situations (b) flow of perfect fluid 
(c) irrotational flows only (d) 2-D incompressible flows. 

(ii) Velocity potential exists for 

(a) all 3-D flow situations (b) flow of perfect fluid 
(c) all irrotational flows (d) steady irrotational flow 

(iii) The continuity equation for a fluid states that 

(a) mass flow rate through a stream tube is constant. 

(b) the derivatives of velocity components exist at every point. 

(c) the velocity is tangential to stream lines. 

(d) Jhe stream function exists for steady flows. 

(iv) V. V = 0 implies that 

(a) dilatation rate for a fluid particle is zero. 

(b) net mass flux from a control volume in any flow situation is zero. 

(c) the fluid is compressible. 

(d) density is a function of time only. 

(v) Momentum theorem is valid only if the fluid is 


(a) incompressible 
(c) inviscid 


(b) in irrotational motion 

(d) irrespective of the above restrictions. 


(vi) Circulation is defined as 

(a) line integral of velocity about any path 

(b) integral of tangential component of velocity about a path 

(c) line integral of velocity about a closed path 

(d) line integral of tangential component of velocity about a closed path. 

(vii) In an irrotational flow, Stokes, theorem implies that circulation is zero 

(a) around two dimensional infinite bodies 

(b) in simply connected regions 

(c) in multiply connected regions 

(d) without any restriction. 

(viii) The curl of a given velocity field indicates 

(a) the rate of increase or decrease of flow at a point 

(b) the rate of twisting of the lines of flow 

(c) the deformation rate 

(d) none of the above 

7.2 Prove that the streamlines «// (r, 6) in a polar coordinate system are orthogonal to 
the velocity potential lines < p(r, 6). 

7.3 The x and y components of velocity in a two-dimensional incompressible flow 
are given by 


u = 3jc + 3y and v = 2x - 3y 


Derive an expression for the stream function. Show that the flow is rotational. 
Calculate the vorticity in the flow field. 


Arts, (-x 2 + 3 xy + (3/2) /, V X V = -1) 
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7.4 A source of volume flow rate 2 m 2 /s is located at origin and another source of 
volume flow rate 4 m 2 /s is located at (3,0). Find out the velocity components at 
(2,2). 

Ans. u — — 0.048 m/s, v — 0.334 m/s) 

7.5 A source of volume flow rate 5 m 2 /s at the origin and a uniform flow of velocity 

4 m/s combine to form a two-dimensional half body. Find out the maximum 
width of the half body. Ans. (1.25 m) 

7.6 A source and a sink of equal volume flow rate 10 m 2 /s are located 2 m apart. If 
a uniform flow of 5 m/s is superimposed, find out the location of the stagna¬ 
tion points. 

Ans. (1.28,0) and (-1.28,0) 

7.7 The discharge of 30 m 2 /s pollutants from a chemical plant into 10 m deep river, 
flowing at 0.3 m/s, can be modelled as a 2-D source across the river depth. It is 
found that the fishes in a certain zone die out whereas those outside the zone are 
unaffected. Find out the extent of this critical zone, if the point of discharge is in 
the midplane of a wide river. 

Ans. (Rankine half body with stagnation point (15.91, 0) and 2v max = 100 m) 

7.8 The Flettner rotor ship (Fig. 7.15) make use of two rotating vertical cylinders. 
Each has a diameter of 3 m and length of 15 m. If they rotate at a speed of 
720 rpm, calculate the magnitude of Magnus force developed by the rotors in a 
breeze of 10 m/s. Assume air density as 1.22 kg/m 3 . 

Ans. (390.083 kN) 

7.9 Find out the strength of a doublet which simulates a physical situation of 2 m 
diameter cylinder in a uniform flow of 15 m/s. 

Ans. (A = 47.124 m 3 /s per metre) 

7.10 Consider a forced vortex rotated at an angular speed ft). Evaluate the circulation 
around any closed path in a forced vortex flow. 

Derive the expression of hydrodynamic pressure as a function of radius for (i) a 
free vortex and (ii) a forced vortex. 

Ans. (Forced vortex Y = 2k r 2 ft); (i) p = - p C 2 /r 2 + C 2 (ii) p = par r 2 /2 + C{) 

7.11 A tornado may be modelled as a circulating flow shown in Fig. 7.24 with 

v r = v z = 0 

Dg = (Or for r < R 


r 



Fig. 7.24 Model of tornado (combination of free and forced vortex) 
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Determine whether the flow pattern is irrotational in either the inner (r < R ) or 
the outer (r > R ) region. Using r momentum equation, determine the pressure 
distribution p (r) in the tornado, assuming p = p^at r —> <*>, Find out the location 
of the minimum pressure. 

Ans. {p =p^ - po? R 2 (1 - r 2 / 2Rr) for r < R; 

p = p„- po? R 4 / 2 r 2 for r > R) 

7.12 A 2 m diameter cylinder is rotating at 1400 rpm in an air stream flowing at 

20 m/s. Calculate the lift and drag forces per unit depth of the cylinder. Assume 
air density as 1.22 kg/m 3 . Ans. (L = 22.476 kN,D = 0) 

7.13 Flow past a rotating cylinder can be simulated by superposition of a doublet, 
a uniform flow and a vortex. The peripheral velocity of the rotating cylinder 
alone is given by v e at r = R (R is the radius of the cylinder). Use the expression 
for the combined velocity potential for the superimposed uniform flow, 
doublet and vortex flow (clockwise rotation) and show that the resultant ve¬ 
locity at any point on the cylinder is given by - 2U 0 sin 0 - v g (at r = R). The 
angle 6 is the angular position of the point of interest. A cylinder rotates at 
360 rpm around its own axis which is perpendicular to the uniform air stream 
(density 1.24 kg/m 3 ) having a velocity of 25 m/s. The cylinder is 2 m in diam¬ 
eter. Find out (a) circulation, F (b) lift per unit length and the (c) position of the 
stagnation points. 

Ans. (236.87 m 2 /s, 7343 N/m, -48.93° and 228.93°) 

7.14 Flow past a rotating cylinder can be simulated by superposition of a doublet, 
a uniform flow and a vortex. The peripheral velocity of the rotating cylinder 
alone is given by v g at r = R ( R is the radius of the cylinder). Use the expression 
for the combined velocity potential for the superimposed uniform flow, 
doublet and vortex flow (clockwise rotation) and show that the resultant ve¬ 
locity at any point on the cylinder is given by — 2U 0 sin 6 - v g (at r = R). The 
angle 0 is the angular position of the point of interest. A cylinder rotates at 
240 rpm around its own axis which is perpendicular to the uniform air stream 
(density 1.24 kg/m 3 ) having a velocity of 20 m/s. The cylinder is 2 m in diam¬ 
eter. Find out (a) circulation, T (b) lift per unit length and the (c) speed of 
rotation of the cylinder, which yields only a single stagnation point. 

Ans. (157.91 m 2 /s, 3916.25 N/m, 382 rpm) 






